Jacobi equations using a variational principle 
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Abstract 

A variational principle is proposed for obtaining the Jacobi equations in sys- 
tems admitting a Lagrangian description. The variational principle gives simul- 
taneously the Lagrange equations of motion and the Jacobi variational equations 
for the system. The approach can be of help in finding constants of motion in the 
Jacobi equations as well as in analysing the stability of the systems and can be 
related to the vertical extension of the Lagrangian formalism. To exemplify two 
of such aspects, we uncover a constant of motion in the Jacobi equations of au- 
tonomous systems and we recover the well-known sufficient conditions of stability 
of two dimensional orbits in classical mechanics. 
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Many classical dynamical systems have a variational formulation, for exam- 
ple, conservative mechanical systems, geodesic flows, classical field theory and even 
geometrical optics. All of them can be described using a Lagrangian function and 
a variant of Hamilton's principle [1]. Given such common factors we call them 
Lagrangian dynamical systems. The variational formulation of such Lagrangian 
systems is not only elegant and compact it also allows uncovering deep connections 
between the dynamical properties and other fields of physics and mathematics [1- 
6] . It is the goal of this short communication to formulate a generalized variational 
principle for Lagrangian dynamical systems, similar to Hamilton's, which is capa- 
ble of producing not only the usual equations of motion but also the dynamical 
equations for deviations between two nearby trajectories: the so-called Jacobi 
variational equations [1,2]. 

We discuss a variational principle for what is known as the equation of geodesic 
deviations in gravitational physics and in Riemannian geometry [2,4,7-8] and give 
a brief outline of the possible applications of its consequences. Let us pinpoint 
that the Jacobi equations can be also regarded as the basic equations for deciding 
questions of dynamical stability, for evaluating the Liapunov spectrum in dynam- 
ical systems [10-11], that they can be shown to occur naturally in the framework 
of the vertical extension of the Lagrangian formalism [3], and, possibly, may be 
used to recast the Jacobi field generated geoodular structure of affinely connected 
manifolds [4]. 

Let us begin formulating the variational principle. Let us assume a iV-degree 
of freedom system that can be described by a Lagrangian L(q, q, £), where q and 
q stand, respectively, for the iV generalized coordinates ((71,(72, • • -q n ) and the iV 
generalized velocities (q\, q%, ■ ■ ■ q n )- Using the Lagrangian of the system, define 
the function 7(q, q, e, e, t) as 

/ • • \ 9L . dL . . 

l{(l,(l,€,€,t) = T —e a + -—e a (1) 
dq a dq a 

here, as in the rest of the article, the summation convention (summing from 1 
to N) is implied for repeated indices. The iV-vectors c = (ei, e 2 , . . . , eAr) and e = 
(ei, €2, . . • , 6n) are to be regarded as describing deviations, and their corresponding 
velocities, from the motion described by q and q. That is, c plays the role of the 
Jacobi field associated with the trajectories of the original system [1,2,4]. Let us 
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notice the important property that 7(q, q, e, c, t) is an explicit function of time 
only when the Lagrangian is time-dependent (i.e. when it is non-autonomous). 

Given the function 7 defined in (1), let us introduce the functional 

E[q(t), £(*)]= / 2 7 (q,q,e,€,t)(i£ (2) 

of the paths joining two given configurations (q\, q\, e\, ei) and (q2, Q2, £2, £2) of 
the varied system between two instants of time t± and £2- The statement of the 
variational principle is now just that 

<5E[q(*),e(*)]=0. (3) 

when the path is varied with the endpoints and the time fixed. The 2N conditions 
guaranteeing that the functional £[(/(£), e(t)] takes an extremum value are the 
associated Euler-Lagrange equations [1] 



= 0, 

(4) 

= 0, a= 1,.. ., AT; 
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or, using the definition (1) in the preceding equations, we obtain the iV Lagrage 
equations of the original system 

d ( dL\ dL „ . , 

= 0, a = l,...,N, (5) 



dt \dq a J dq a 
plus the iV equations: 

M ab e b + Cabh + K ab e b = 0, a = l,...,N, (6) 

for the deviation, e, between two nearby trajectories. The N x N matrices M, C 
and K, are defined by 



M ab 
K ab 



d 2 L 

dq a dq b 
d_ f d 2 L 
dt \q a qb 



Cab = 

d 2 L 
dq a dq b 



\- 


( 92L ) 


dt 


\QaQbJ 



a, b 



+ 



d 2 L 



dq a dq b 
1,...,N. 



d 2 L ' 
dq b dq a _ 



(7) 
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Equations (6) are the Jacobi variational equations for the original system [1,2,4]. 

If the Lagrangian of the system is time-independent, the system has a well- 
known constant of motion 



in such time-independent case the Jacobi equations also admit — hence the impor- 
tance of the property of 7 mentioned above — an analogous constant, namely 

h = -^qb + ^e 6 ; (9) 
oq h oe b 

using definition (1), the constant h can be cast in the form 

. dH . OH 

h =w^w^ (10) 

Equation (9) [or (10)] is an important conclusion from the variational formulation. 

Notice that the variational equations can be written in first-order form pro- 
vided that 

HaS^h <"> 

i.e. that the matrix M is invertible; this is also the condition for the existence 
of a Hamiltonian description of the system [1]. Granted such condition, the 2N 
variational equations can be written as 

where the 2iV-vector x and the 2N x 2N matrix J7", are defined respectively by 
x= (i)' ^ = \ A/r-iv ns-ir> ) > ( 13 ) 



JJ' ° \M- X K M~ X C f 
where the and the 1 stand, respectively, for iV x N zero and unit matrices. 

We can now, for example, use the solutions to equations (12) [or to (5) and 
(6)] to evaluate the iV Liapunov exponents, A a , in the standard way [8,11] 

A a = hm \ log H^fI, a=l,...,N. (14) 
t-00 t ||x a (0)|| 
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where ||v|| is any norm of the iV-vector v and x(0) is just the initial condition 
for a perturbation, in one of the N appropriate directions, of the orbit q(i) under 
analysis. 

We emphasize that the variational formulation of the Jacobi equations can be 
of help for discovering constants of motion in the variational equations of partic- 
ular systems (using their symmetries, for example), like the one we derived here 
[equation (9)] for the autonomous case. The close relationship of the variational 
principle to the Lagrangian through definition (1) reflects the central relation- 
ship between L(q, q, t) and the properties of the motion. The Lagrangian is thus, 
paraphrasing an apt description, the true gene of the motion [12]. 

The results derived from our variational principle have some bearings on mat- 
ters of stability. For example, using our results, we can analyse the stability of 
particle orbits in two dimensions recovering well-stablished results in the process. 

To illustrate the previous assertion, let us analyse the motion of a particle with 
unit mass under a particular time independent potential, U (q) , in two dimensions 
[9,13]. Select a particular orbit O as a reference, the generalized coordinates can 
now be chosen as the distance, z, from the reference orbit to the particle and the 
arc length, s, from an arbitrary origin on O to the point, also on the reference 
orbit O, from which z is measured. The Lagrangian of the system and the function 
7 in these coordinates are (in obvious notation) 



7 ( e s j ; ) 




) 



= ze z + se s 1 + - + -2 
V PJ LP 

zs 2 dp dU 




dz 



dU 



(16) 




where p is the radius of curvature of the reference orbit at s. The system is 
autonomous, hence the quantity 



is a constant of the motion. Combining (16) and (17), using equations (5) and 
(6), and choosing O as the orbit under analysis, we can show that the equation 
governing the stability of that orbit is 



since the e s equation concerns itself only with deviations tangent to the original 
orbit. The stability of the two-dimensional orbits is easily established in the h = 
case — which simply corresponds to analysing varied trajectories in which the 
energy does not change respect to its value in the nonvaried orbit O. In this 
homogeneous well known result [14] asserts that the above equation has 

oscillatory solutions (hence, the analysed orbit is stable) if the quantity between 
parenthesis in (18) — called the coefficient of stability — is positive definite at every 
point on the orbit [13]. 

For the relationship to the stability of periodic orbits see, for example, [13, 15]. 
Furthermore, we think the variational principle (3) can be useful for explaining in 
non-standard way some results concerning the relationship between singularities 
in the exponential map and the corresponding Jacobi fields in geodesic flows on 
manifolds [2] with its natural geoodular structures [4], and can have relevance 
in the study of both Lagrangian and time-dependent Hamiltonian mechanics and 
field theories in the context of their vertical extensions [3] . 
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